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ABSTRACT. A significant fraction of mass in the universe is believed to 
be in the form of dark matter (DM). Due to gravitational instability, the 
DM collapses hierarchically into DM halos. In this contribution we present a 
study of the formation and evolution of such DM halos in a COBE-normalized 
spatially flat ACDM model (Qo = 1— Qj^ = 0.3; h = 0.7) using high-resolution 
iV-body simulations. The novelty of this study is use of the newly developed 
halo-finding algorithms to study the evolution of both isolated and satellite 
(located inside virial radii of larger group- and cluster-size systems) halos. 
The force and mass resolution required for a simulated halo to survive in the 
high-density environments typical of groups and clusters is high: ~ 1 — 3 kpc 
and ~ 1O 9 M0, respectively. We use the high-resolution Adaptive Refinement 
Tree (ART) ./V-body code to follow the evolution of 256 3 dark matter particles 
with dynamic range in spatial resolution of 32, 000 in a box of 60/i _1 Mpc. 
We show that the correlation function of these halos is anti- biased with respect 
to the dark matter correlation function and is high and steeper than the 
correlation function of the isolated virialized objects. The correlation function 
evolves only mildly between z = 3 and 2 = 1. The mass evolution of isolated 
virialized objects determined from the simulation is in good agreement with 
prediction of semi-analytical models. The differences exist, however, if we 
include satellite halos in the halo catalogs. 



1 Introduction 

It is generally believed that dark 
matter (DM) constitutes a large 
fraction of the mass in the Uni- 
verse. Therefore, it significantly af- 
fects both the process of galaxy for- 
mation and the large-scale distribu- 
tion of galaxies. The most convinc- 
ing observational evidence for sub- 
stantial amounts of dark matter even 
in the very inner regions of galax- 
ies comes from recent HI studies 
of dwarf and low surface brightness 
galaxies. The gravitational domina- 
tion of DM on the scale of galaxy 
virial radius implies that collision- 
less simulations can be used to study 
formation of the DM component of 
galaxies. 

A well known problem of dissi- 
pationless simulations is overmerg- 
ing, i.e. the lack of substructure in 
virialized objects that could be as- 
sociated with galaxy locations. This 
effect is due mainly to insufficient 
force and mass resolution (Moore et 



al. 1996; Klypin et al. 1998, here- 
after KGKK) . Recently, the dynamic 
range of the iV-body simulations has 
become sufficiently high to overcome 
this problem. With sufficient reso- 
lution the DM halos survive and 
can be identified in high-density re- 
gions even in collisionless simulations 
(KGKK; Ghigna et al. 1998). This al- 
lows us to study halo dynamics and 
the physical effects of tidal stripping 
and dynamical friction in group- and 
cluster-size objects. 

In this contribution, we discuss 
the evolution of DM halos in a spa- 
tially flat cosmological model domi- 
nated by cold dark matter and non- 
zero cosmological constant (ACDM). 
Specifically, we study the differences 
between isolated and satellite halos 
in their mass evolution and spatial 
distribution. Scenarios with a A-term 
have become very successful in de- 
scribing most aspects of cosmological 
structure formation. We have cho- 
sen the following parameters in or- 
der to reconcile the model with both 
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Figure 1. The correlation function of halo s with circular velocities v c i rc > 120 
km/s (solid line) and u c ; rc > 150 km/s (dashed line) in comparison to the 
correlation function of dark matter particles (dotted line). The dashed-dotted 
lines represent the correlation function of objects found by the friends-of- 
friends algorithm with the linking radius corresponding to the virial overden- 
sity (dash-dotted for objects with N > 100 particles and dash- triple- dotted 
for N > 20). 




the COBE measurements and clus- 
ter abundance: f2o = 1 — Ha = 0.3; 
h = 0.7; erg = 1.0. The age of the uni- 
verse in this model is ~ 13.5 Gyrs. 

2 Numerical Simulation 

In order to study the properties of 
halos in a cosmological environment, 
the simulation box should be suffi- 
ciently large. On the other hand, to 
assure that halos do survive also in 
clusters, the force resolution should 
be ~ 1 — 3fe~ 1 kpc and the mass 
resolution should be < 1O 9 /i _1 M0 
(Moore et al. 1996; KGKK). The 
Adaptive Refinement Tree (ART) 
Af-body code (Kravtsov, Klypin & 
Khokhlov 1997) reaches a formal dy- 
namical range of 32, 000 in high den- 
sity regions, which for the 60 h^ 1 
Mpc box corresponds to the required 
necessary force resolution. In the 60 
h^ 1 Mpc box with 256 3 particles, 
each particles has a mass of 1.1 X 
1O 9 /i- 1 M . 

Identification of halos in dense 
environments and reconstruction of 
their evolution is a challenge. Most 
widely used halo-finding algorithms, 
the friends-of-friends (FOF) and the 
spherical overdensity, both discard 
"halos inside halos", i.e. satellite ha- 
los located within virial radius of 
larger halos. The distribution of ha- 
los identified in this way, cannot be 



compared easily to the distribution 
of galaxies, because the latter are 
found within larger systems. In order 
to cure this, we have developed two 
related algorithms, which we have 
called hierarchical friends-of-friends 
(HFOF) and bound density maxima 
(BDM) algorithms (KGKK). These 
algorithms are complementary. Both 
find essentially the same halos. The 
advantage of the HFOF algorithm 
is that it can handle halos of ar- 
bitrary (not only spherically sym- 
metric) shape. The advantage of the 
BDM algorithm is that it separates 
background unbound particles from 
the particles gravitationally bound 
to the halo, and thus allows a better 
determination of physical properties 
of halos. 

Since the algorithms work on a 
snapshot of the particle distribu- 
tion, they tend to identify also small 
fake "halos" consisting of only a 
few unbound particles, clumped to- 
gether by chance at the analyzed mo- 
ment. We deal with this problem by 
both checking whether the identified 
clump is gravitationally bound and 
by following the merging history of 
halos. Halos that do not have a pro- 
genitor at a previous moment are 
discarded. For other halos we find 
the direct progenitor, i.e. a halo at 
a previous moment that contains the 
maximum number of particles of this 
halo. We use the chain of progenitors 
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identified in this way to reconstruct 
the mass evolution of a given halo 
back in time, down to the epoch of 
its first detection in the simulation. 

About 10,000 halos (with maxi- 
mum circular velocity > 120 km/s, 
cf. Kravtsov et al. in this volume) can 
be identified at z = in the analyzed 
60 h~ 1 Mpc simulation box. Several 
hundreds of these halos are located 
in groups and clusters. This allows us 
to carry out for the first time a sta- 
tistical comparison study of the clus- 
tering and mass evolution of isolated 
and satellite galaxy-size halos in an 
hierarchical cosmological model. 

3 Results 

We have identified halos using dif- 
ferent input numerical parameters 
(size and number of particles) for 
the halo-finding algorithms. We find 
that for distances > 100/i — 1 kpc the 
resulting halo-halo 2-point correla- 
tion function (CF) does not depend 
on these assumptions, as long as 
the number of particles in identi- 
fied halo is > 30. Note that the halo 
radius does not set limits on the 
inter-halo separation because halos 
are allowed to overlap. In Fig. 1 
we present the halo-halo correlation 
function for halos with a maximum 
radius of lOO/i^ 1 kpc and more than 
28 particles (3 X lO lo fe _1 M ). There 
is a statistically significant anti-bias 
on scales less than the correlation 
radius ro- The slope of the CF is 
fa —1.8 at scales fa 1 — 15fe _1 Mpc; 
at smaller scales the correlation func- 
tion flattens slightly. In agreement 
with Colin et al. (1998) and results 
of many other studies, we find almost 
no time evolution of the CF at z < 3. 
The correlation function of DM, on 
the other hand, evolves rapidly which 
results in the strong evolution of bias 
and transition from b > 1 to b < 1 
at small scales at z fa 1 (sec Fig. 2, 
see also Kravtsov et al. on the bias 
evolution in this volume). 

The correlation function for the 
halos identified using our algorithms 
is compared in Fig.l to the CF of 
halos identified in a standard way 
using the FOF algorithm with link- 
ing radius corresponding to the virial 
overdensity. As we noted above, by 
definition the satellite halos are not 
present in the catalogs so gener- 
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Figure 2. The correlation function of 
halos with circular velocities v c i rc > 
120 km/s at z = (solid line), 
z = 1 (dash-dotted line), and z = 3 
(dashed line) in comparison to the 
correlation function of dark matter 
particles at the same redshifts (dot- 
ted lines; amplitude increases with 
decreasing redshift). 

ated. The comparison shows clearly 
that there are significant differences 
in both the slope and the ampli- 
tude at small (<7/i Mpc) scales 
between these correlation functions. 
The higher amplitude of the CF de- 
termined using our halo-finding al- 
gorithms is explained by the larger 
number of small-separation halo 
pairs formed by satellite halos in 
groups and clusters missing in the 
FOF catalog. Note also that the 
CF of FOF halos docs not extend 
down to 200h _1 kpc due to the self- 
exclusion of halos at separations 
smaller than halo virial radii. 

The mass of an object found by 
the HFOF algorithm at virial over- 
density can be defined as the sum of 
linked particle masses. For all of the 
HFOF objects we identify the main 
progenitors at all epochs down to the 
halo formation time. To study the 
mass evolution due to merging and 
accretion we have divided these ob- 
jects into four mass bins at z = 0: 
M > 5 x 10 12 M Q (bin 1), 5 x 
1O 12 M > M > 5 x 10 n M Q (bin 
2), 5 x lO n M > M > 5 x 10 10 M Q 
(bin 3), and M < 5 X 10 10 M Q (bin 
4). Average mass evolution of halos 
in these bins normalized to the mass 
at z = is shown in Fig. 3. We 
find a good agreement with the semi- 
analytical predictions (Lacey & Cole 
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Figure 3. Mass evolution of halos 
identified by the HFOF algorithm at 
virial overdensity. The dotted line 
is for average halo mass of 1.4 X 
1O 13 M (10 halos), the dashed for 
1.1 X 10 12 M Q (20 halos), the dash- 
dotted for 7.7 x 10 10 M Q (417 ha- 
los) and the dash-triple- dotted line 
is for 2.1 X 10 10 M Q (82 halos). The 
solid lines show the predictions of 
semi- analytical model (kindly pro- 
vided by Claudio Firmani) for 10 13 , 
10 12 , and lO n M , from the bottom 
to the top. 

1993) for the evolution of these iso- 
lated halos. 

Unfortunately, there is no simple 
and straightforward way to assign a 
mass for all halos identified in the 
simulation. Unlike the isolated halos 
identified by HFOF at virial overden- 
sity, the satellite halos, although sur- 
viving, are subject to the tidal strip- 
ping which reduces their mass. They 
are limited therefore by tidal, rather 
then virial, radius. To assign masses 
to the halos we proceed as follows. 
The isolated halos are assigned the 
mass inside the virial radius or radius 
of lOO/i -1 kpc, whichever is smaller. 
The satellite halos are assigned the 
total mass of gravitationally bound 
particles within their tidal radius (or, 
again, within 100/i -1 kpc, whichever 
is smaller). The tidal radius is deter- 
mined as the radius at which the den- 
sity profile of a halo flattens (stops 
decreasing). 

We now construct the complete 
mass evolution histories for the set 
of all halos with the masses assigned 
as described above. We have divided 
these halos into five groups with 



masses M > 1O 13 M , 1O 13 M > 
M > 5 x 1O 12 M , 5 x 10 12 M Q > 
M > 10 12 M Q , 10 12 M Q > M > 5x 
10 n M Q , and 5x 10 n M Q > M . We 
defined a subset of 3674 halos, mass 
of which increases (with allowance 
for small statistical fluctuations) at 
all epochs. As before, the mass of 
these objects is normalized to their 
final mass at z = 0. The mass evolu- 
tion of these halos is shown in Fig. 4 
(solid lines). The overall evolution is 
similar to the mass evolution of iso- 
lated halos described above (Fig. 3). 
Note, however, that while the mass 
evolution tracks are curved in Fig. 3, 
the mass evolution of the sample that 
includes satellites can be better rep- 
resented by the straight lines in these 
log-log plots. This difference is due 
to the different halo selection proce- 
dure and to the different assignment 
of mass to the selected halos. 

In the two lowest mass ranges we 
also find an additional subset of 2650 
halos, whose masses decrease after 
z = 1 (dashed and dot-dashed curves 
in Fig. 4). Their mass increases at 
high redshifts, reaches a clear max- 
imum and decreases thereafter. The 
mass of these halos grows first due 
to accretion of surrounding material 
and of smaller halos. At some point, 
however, these halos are being ac- 
creted by more massive halos and 
they start to loose mass due to the 
tidal stripping and interaction with 
other satellite halos. At z = 1 these 
halos arc distributed similarly to the 
the rest. At z = 0, however, they 
are clustered more strongly than the 
overall halo population. 

This is illustrated in Fig. 5, which 
shows the correlation function for 
the two subsets of halos: always in- 
creasing mass and decreasing mass 
at z < 1. The correlation functions 
of the former has a lower amplitude 
and is not as steep as the corre- 
lation function of the latter. Note 
that the CF of the halos with the 
ever-increasing mass is anti-biased at 
scales < lOh^ 1 Mpc, the CF of the 
halos that loose mass is actually pos- 
itively biased. This reflects the fact 
that the loosing mass halos arc found 
within massive systems such as mas- 
sive galaxies, groups, and clusters, 
and are therefore strongly clustered. 

One might speculate that this dif- 
ference in the correlation functions 
may serve as a possible cxplana- 
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Figure 4. Mass evolution of ha- 
los. The solid lines are for average 
masses of (from the bottom to the 
top) 1.2 X 10 13 M Q (14 halos), 6.6 X 
10 12 M Q (34 halos), 1.9 x 10 12 M Q 
(442 halos), 7.0 X IO^Mq (534 ha- 
los), and 2.4 X 10 11 Mq (2650 ha- 
los). The dot- dashed (average mass 
o/6.9x IO^Mq) and the dashed (av- 
erage mass of 2.0 X 10 11 Mq) lines 
show the mass evolution of a sub- 
set of halos which loose mass between 
z = 1 and the present due to the tidal 
stripping in groups and clusters. 



tion for the color segregation of the 
correlation amplitude that has been 
recently observed (Carlberg et al. 
1998). In fact, one could expect that 
the galaxies hosting halos which un- 
dergo different mass evolution also 
show different properties, and col- 
ors in particular. Further studies are 
necessary to test whether this simple 
picture can really explain the obser- 
vations. 
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Figure 5. The correlation func- 
tion of halos with circular velocities 
Vcirc > 120 km/s. The solid line cor- 
responds to halos the mass of which 
always increases (solid lines in Fig. 
4), the dashed line corresponds to 
halos which loose mass during evolu- 
tion (dashed and dashed- dotted lines 
in Fig. 4)- The correlation function 
of dark matter particles is shown by 
the dotted line. 
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